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5-H ' Abstract 

H ' Two channels of the sd exchange interaction are considered in magnetic 

^ ; , junctions. The first channel describes the interaction of transversal spins 

. with the lattice magnetization. The second one describes the interaction of 

C ' longitudinal spins with magnetization. We show the longitudinal channel 

I , leads to a number of significant effects: 1) drastic lowering of the current 

' instability threshold down to three (or even more) orders of magnitude; 

^ ' 2) creation of large enough distortion of equilibrium due to current driven 

O . spin injection leading to inversion of energy spin subband populations and 

, ^ , ' laser-like instability in THz frequency range at room temperature. 

I External magnetic field may tend to lower additionally the instabil- 

T-H ■ ity threshold due to the proximity effect of purely magnetic reorientation 

^ ' phase transition. This effect demonstrates the new properties: the gi- 

, ant magnetoresistance (GMR) becomes strongly current dependent and 

• the exchange switching becomes of very low threshold. We derived some 

nI ' matching condition that should be satisfied to achieve high spin injection 

, level. Some characteristic quantities were appeared in the condition. We 

' investigated also the junctions having variable lateral dimensions of the 

T-H I layers, for example, a ferromagnetic rod contacting with a very thin fer- 

. romagnetic film. Large enhancement of the current density may appear 

' near the contact region leading to the spin injection luminescence. 

> : 
•1—1 , 

^ ■ 1 Introduction 

5t 1 The last years a spin dependent transport in ferromagnetic films and junc- 
tions becomes of growing interest. The first question arising in this field 
is the nature of the interaction between conduction electrons and lat- 
tice magnetization. The simplest approach proposed by Vonsovskii [T] is 
known as the sd exchange model. As it will be shown below, two dif- 
ferent channels of the sd interaction may be detailed. The first channel 
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describes the transversal spin transformation leading to the spin transfer 
torque, while the second channel describes the longitudinal spin transfor- 
mation leading to distortion of spin subband populations. We sketch here 
the both channels in the frame of a unique theory and show the main 
peculiarities (for more detail see Refs. [H [3] S) [5]). 

Historically, the large significance had the so called " Giant Magnetore- 
sistance" (GMR) effect when reorientation phase transition was observed 
in ferromagnetic junction at some critical external magnetic field He- In 
modern experiments, the value of this effect may be large enough, up to 
tens percentages. This effect is good developed now and the Nobel Prize 
was awarded in 2007 to physicists A. Pert and P. Griinberg for discovery 
of the effect [IE]. 

The next researches in the field concerned the current flowing in the 
junction. As it was predicted [SlIU, magnetization instability may arise 
at current density exceeding some threshold value, j > jth ~ lO'^-lO* 
A/cm^. The transversal channel of the sd exchange is responsible for 
the instability. Many experimental confirmations of the instability were 
obtained starting from the first one [10]. A principal problem remains, 
however, namely, a relatively large current threshold. It would be very 
interesting to estimate the threshold for the second longitudinal channel 
of sd exchange interaction. We try to answer this question below. 

2 The structure investigated. Mechanisms 
of exchange switching 

We take as a starting point the simplest plane structure containing two 
ferromagnetic layers shown in Fig. 1. The layer 1 has pinned lattice 
magnetization Mi, while layer 2 has free lattice magnetization M; by a 
convention, arrows denote both magnetizations. Conduction electrons, of 
course, have free spins everywhere. The electron current density j/e flows 
perpendicular to the layers, e is the electron charge. 

Two switching mechanisms may be seen from the Fig. 2 (see also [S]). 
The first one arises due to noncoUinearity of the vectors Mi and M and 
the loss of the transversal spin components during their moving in the 
layer 2. The lost components transfer from the mobile electrons to the 
lattice [8l[9], which may excite the switching. Then the electron spins be- 
come completely coUinear with M, but remain nonequilibrium ones. The 
second mechanism arises at this stage. Equilibrium distribution (among 
the spin energy subbands) should be restored and processes go, which may 
lead to M instability and switching [111 112) . 

3 Equations 

We intentionally exclude from consideration some zone between the layers 
1 and 2 where a quantum nonuniform problem should be solved to describe 
the conduction electrons moving between the layers. Instead, we tried to 
derive some flux boundary conditions which allow considering the junction 
processes without the detail inside the zone (all the theory was presented 
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Figure 1: Plane scheme of the magnetic junction. There exists a very thin 
spacer between the layers 1 and 2; H is an external magnetic field; the layer 3 
is nonmagnetic. 




m ^0,m^0 m =0 



Figure 2: Scheme of processes in the layers of the junction. 
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in Refs. [2][3l|4])). The approaches of some other authors are presented 
in Refs. [HI [SI HI]. 

The standard sd Hamiltonian is 

Usd = —cti / mi (x')Mi da;' — Q / m(x')M(a;') da;', (1) 
J -L Jo 

where q ~ qi ~ 10"* is the dimensionless sd exchange parameter, Li ^ 
L. The dynamics may be described by the following equations: 

1) the continuity equation for mobile electrons 

9m dJ . , m — m „ 

_ + _+^«[„,xM] + -— =0, (2) 

where the spin flux density for mobile electrons J is 

J = ^(jt-iJM, (3) 

and j_i being the partial currents, M — M/|M| is the unit vector, r is 
the spin relaxation time, pi_g is the Bohr magneton, 7 is the gyromagnetic 
ratio, m is the equilibrium value; 

2) the Landau-Lifshitz-Gilbert equation for M 



— + 7 [M X Heff J - K 



(4) 



where k is the Gilbert dissipation constant (typically k 
Effective field in Eq. ((U is 

Heft =H + H„+Hd+H,d + ^-^, (5) 

where H is the external magnetic field. Ha is the anisotropy field, Hd is 
the demagnetization field and H^^ is the sd exchange field. Here all the 

SU 

fields are defined by external conditions, except Usd = which is a 

functional derivative and should be calculated from Eq. ([ij. We suppose 
very small spin relaxation time r ~ 3 x 10~^^ s, so that we have ujt ^ 1 
condition for the characteristic precession frequency u. Along with it, 
specific exchange frequency yaM is large enough, namely, ujsd ~ 10^ s~ , 
and therefore cjsdT ~ 10^ 2> 1 . Based on the assumptions mentioned, we 
may solve Eq. and substitute the solution into Eq. U]). It was exactly 
performed in Refs. [3] [H [5] with the result for the nonequilibrium part of 
the exchange field Had 

AUsd^hsdMil5{x^0), (6) 

where field h^d is a function of the conduction electron parameters de- 
pending on the current direction (forward or backward), I is the electron 
spin diffusion length in layer 2, and the 5-function shows the field © 
is localized near the a; = boundary of the layer (see the derivation in 
Refs. [31 [H [5]). Thus, we present the form of equation ([6]) and now it is 
necessary to formulate the boundary conditions to solve the problem. 
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We formulate further the conditions of magnetic flux continuity fol- 
lowing the results of Refs. I31|31[S]. 

The following magnetization flux densities exist in our problem [3]: 

1) Free electron spin current density J (see ([3])). It is a longitudinal 
flux because J ~ M. 

r - dM 

2) Lattice magnetization flux density Jm(x) = 7AM M x — — 

ox 

which is, obviously, transversal. 

3) Sd exchange flux density Jsd{x) = jhsJ [m(0) x Mij e(x - 0), 

which is transversal, 6(x) = 1, a; > and ^(a::) = 0, a; < 0. The general 
continuity condition will be 

J(+0) - J(-0) + Jm(+0) - Jm(-O) + J,d(+0) - J,d(-0) = 0. (7) 

The situation may be simplified for pinned layer 1 when Ja/(— 0) = 
Jsd(— 0) — . Moreover, let us consider separately the projections for 
the forward and backward currents [51[S]. Then we have for the forward 
current {j/e > 0) 

[m(+0) X [j(-0) X M(+0)]] = Jm(+0) + J.d(+0), (8) 

and for the backward current (j/e < 0) 

[Mi X [j(+0) X Ml]] = -Ja/(+0) - J.d(+0). (9) 

In the other boundary of layer 2, that is at a; = L, we have 

Jm(L-0)=0. (10) 



4 Instability thresholds 

We solve now Eqs. ([2]) and Q with boundary conditions ((8))- (|10p . Initial 
magnetization is taken as Mz = ±1, that is directed along z axis. We 
search the small harmonic fluctuations AMx, AMy ~ exp{—iLjt) and flnd 
dispersion relations (see [3l|4l[5]). Then we obtain the following estimation 
of the instability threshold currents. Threshold current density jth for the 
transversal channel is 



Jth 
e 



2.,M^lX.f^^^ (11) 



where A = L/l, Qi is the current polarization degree in layer 1, v* = 
Z1/Z3 + \Z\IZ2, Zi (i = 1, 2, 3) is the spin resistance (see [5]). Estima- 
tions using Eq. (fTTj) typically give jth ~ 6 x 10^ x ^1 -I- A/cm^. 
Threshold current density for the longitudinal channel is 

,M _ m ^+|)(l + ^g.), (12) 
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where we may substitute the typical values Ha ~ 100 Oe, A « 0.1 ^ 
Z2/Z3, I ~ 10~^ cm. Then at iJ" = we obtain the threshold jth ~ 
2 X 10'' A/cm^ which is three orders of magnitude low than for transversal 
channel. Further radical reduction of the threshold may be achieved in 
a magnetic field HMz — i- —Ha + 0. The cause of the reduction is the 
proximity to the reorientation phase transition in the field H. 

5 Rod-to-film cylindrical structure 

Up to now, we considered the simplest planar ferromagnetic structures. 
However, to obtain the highest spin injection level we should check other 
configurations also. Further we will focus our attention on a cylindrical 
structure of the rod-to-film type, the scheme of which is shown in Fig. 3. 
This scheme was proposed in Ref. [15|. The thickness h of the film 1 is 
small in comparison with the radius R of the rod 3. Due to continuity 
of the current we may get R/2h ~ 500 times enhancement of the current 
density near the edge of the rod when the current flows from the rod to the 
film. Then the current density may reach 10^ A/cm^ or even more. It is, 
apparently, enough to have very high spin injection level and the inversion 
of spin population. The possibility of such an inversion was discussed in 
a number of papers [161 1171 1181 119] . We consider below some calculations 
and experimental results about spin-injection and THz luminescence in 
the structures. 

The distribution of electron spins in the structure under electron cur- 
rent flowing in the rod — >■ film direction is calculated by means of conti- 
nuity equation (see [15| ) 

y^p~^^-^-o, (13) 

where P = (n-f — ni)/n is the degree of spin polarization and rif, ni are 
the populations of the lower and upper spin energy subbands, n = n-|- + nj_, 
P is the equilibrium polarization, jo = enl/r is a characteristic current 
density. 

We solve Eq. (|13p analytically in cylindrical coordinates using the con- 
ditions of spin flux continuity at the layer boundaries. The polarization 
tends to its equilibrium value P, when we are moving apart from the rod 
edge. 

The calculated P{r)/P distribution is shown in Fig. [4] Curves 1-4 
correspond to rising of the spin injection level Q\/P: 0(1), 1(2), 2(3), 
5(4). We see an inversion of spin polarization, P{R) < appears for large 
enough value of Qi/P. 

The calculated dependence of the relative spin polarization P{R)/P 
on the dimensionless current density j{R)/jD is shown in Fig. O As it is 
seen, the inversion of spin population {P{R) < 0) appears also. However, 
for a nonmagnetic rod (curve 1, Qi/P = 0) the inversion is absent. It 
means the inversion is the consequence of spin injection by the current. 
The negative polarization rises in magnitude with current growing, that 
is, with the injection level. 
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10 nm 



Figure 3: Structure scheme: 1-ferromagnetic film, 2-substrate, 3-ferromagnetic 
rod, 4- nonmagnetic conductor, 5-fluoroplastic plate, U voltage. 
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Figure 4: The calculated P{r)/P distribution is shown in the figure. The fol- 
lowing parameters are taken: R/l = 20, j{R)/jD = 1 and various values of 
Qi/P ratio. We see an inversion of spin polarization, P{r) < 0, appears for 
large enough value of Qi/P- 
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Figure 5: Spin polarization at the boundary between the rod and the film as a 
function of the (dimensionless) current density j{R)/jD a.t R/l — 20 and various 
values of Qi/P ratio. 

6 Experimental observation of spin- inject ion 
luminescence 

We measured the luminescence of the experimental structure (Fig. 3) by 
means of a Golay cell and the metallic filter to cut off frequencies below 
approximately f THz. 

Measurements of luminescence intensity were carried out with the 
pulse current flowing in forward and backward directions. Pulses may 
be of different pulse period to pulse duration ratios (PPPDR) to have 
small and variable heating. 

The results are presenting in Figs. |6] and [T] We see the measured 
intensity depends on the current direction. These observations cannot be 
explained by any thermomagnetic effects, such as Peltier or Ettingshausen 
effects [20]. For metals, the latter effects may be estimated as a fraction 
of a degree, while we have heating up to 10-^15 degrees. Moreover, the 
discussed dependence on current direction disappeared immediately after 
we replaced the steel rod by the nonmagnetic copper one. We conclude 
therefore, the direction of the current flowing influences due to magnetic 
properties of the rod and represents non-thermal action of the current. 

Based on the observations, we suggested the main role of magnetic 
properties of the rod and non-thermal action of the current. The lumines- 
cence really contains thermal and non-thermal contributions. However, 
the thermal contribution decreases with PPPDR increasing, while the 



9 



1.0 - 




Figure 6: Measured intensity / of the structure radiation intensity as a function 
of time for the forward (1) and backward (2) current directions, PPPDR being 
5. Initially / rises with time t due to the current heating and spin injection. 
After the current is switched off, the radiation drops. 
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Figure 7: Measured intensity / for the forward (1) and backward (2) current 
directions, PPPDR being 20. Thermal effect is greatly suppressed, but the 
difference between forward and backward current effects increases. 



non-thermal contribution remains stable. That is the main cause of the 
splitting curves 1, 2 in Figs. [6] and [7] 

We may use now the results of spin-injection calculations derived re- 
cently in p3]. According to Eq. (20) from the Ref. :15^, we represent the 
highest (in magnitude) negative value of the nonequilibrium spin polar- 
ization achieved in film near the boundary of the rod in the form 



1AP|= |Qi(Mi.M(i?))-p]^'(^) ^-(^Z^) 



(14) 



JD K.+iim 

where K,j is the modified Bessel function of the second kind with index 

^4?^. (15) 

2 I jD 

The most significant consequence of the formulae (|14p and (I15[l is the fact 
that the nonequilibrium polarization AP depends on the current both di- 
rectly and via the index v, being nonsymmetrical with respect to changing 
the current sign, j — > —j. Therefore the spin injection contributes to AP, 
and the contribution is different for forward and backward directions of 
the current. The absolute difference between the contributions, accord- 
ing to the formulae, have no small parameters and, in principle, may be 
sufficient to explain the splitting 1, 2 curves in Figs. [S] and [T] 
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7 Summary 



The conduction electrons that participate in polarized current {s elec- 
trons) interact with the lattice magnetization (rf electrons) in a ferro- 
magnetic junction via two channels: (i) via the transfer of the transverse 
spins (perpendicular to the lattice magnetization) to the lattice, and (ii) 
via the transfer of the longitudinal spins parallel to the magnetization to 
the spin energy subbands. The latter can be considered as a change in the 
population of spin energy subbands, i.e., the injection of nonequilibrium 
spins. This injection leads to the creation of a nonequilibrium sd exchange 
effective field, which, in turn, affects the dynamics of the system, in par- 
ticular: 1) the lowering of magnetic exchange instability threshold, and 2) 
creation of the inversion subband population and negative effective spin 
temperature. 

We should provide some specific relations between the spin resistances 
of the layers Zi, where i — 1, 2, 3 labels the layers. In particular, Z2 <^ 
Zi, Zz condition leads to reduction in the threshold. The estimates for 
some typical samples show the threshold can be lowered by orders of 
magnitude, for example, from ~ 6 x 10^ to 2 x 10^ A/cm^. The minimum 
thresholds always correspond to the predominance of the spin-injection 
channel of the sd exchange interaction. 

An external magnetic field H which is near the critical value He for a 
reorientation phase transition {H — He — 0) can lead also to radical lower- 
ing of the exchange current threshold. The external magnetic field, being 
near the phase transition point and acting together with the exchange 
field, helps the exchange switching. We investigated also the junctions 
having variable lateral dimensions of the layers, the so called rod-to-film 
structures. Very high current density and spin-injection level may be 
achieved in the structures. Two interesting facts have been observed in 
our measurements: 1) the presence of non-thermal contributions to TIfz 
luminescence from the system in study, and 2) the difference between the 
radiation intensities under the forward and backward current directions. 
As it was shown (see Eqs. p4|) and H15|) '). the spin injection in the junction 
depends substantially on the current direction. Therefore, the facts men- 
tioned may be due to the radiation created by the nonequilibrium spins 
injected near the rod. 
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